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Abstract 

We calculate the total cross section for the production of photons from the scattering of axions 
by a strong inhomogeneous magnetic field in the form of a 2D (5-function, a cylindrical step function 
and a 2D Gaussian distribution, which can be approximately produced by a solenoidal current. 
The theoretical result is used to estimate the axion-photon conversion probability which could be 
expected in a reasonable experimental situation. The calculated conversion probabilities for QCD 
inspired axions are bigger by a factor of 2.67 (for the cylindrical step function case) than those 
derived by applying the celebrated ID calculation of the (inverse) coherent Primakoff effect. We 
also consider scattering at a resonance Eaxion ^ "^axion: which corresponds to the scattering from 
a (5- function and gives the most enhanced results. Finally, we analyze the results of this work in 
the astrophysical extension to suggest a way in which they may be directed to a solution to some 
basic solar physics problems and, in particular, the coronal heating problem. 
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I. INTRODUCTION 



The possible existence of a light pseudoscalar particle is a very interesting possibility. For 
example, the axion , which was introduced in order to solve the strong CP problem in 

QCD, has since then also been postulated as a candidate for the dark matter in the universe. 
A great number of ideas and experiments for the direct detection of this particle have been 
proposed in the past jj], jsl. For example, It was recognized by Sikivie that axion detection 
exploiting axion to photon conversion in a magnetic eld was a possibility (6|. 

Related to that, in a series of recent publications by one of us fz^, it was shown that an 
axion-photon system displays a continuous axion-photon duality symmetry when an external 
magnetic field is present and when the axion mass is neglected. This allows one to analyze 
the behavior of axions and photons in external magnetic fields in terms of an axion-photon 
complex field. For exarnple, the deflection of light from magnetars has been recently studied 
using these techniques jsj. It is important to note here that the same duality symmetry exists 
also when considering massive photons, under the condition = rria, that is the photon and 
the axion masses are equal. These conditions can be achieved when conducting experiments 
where the axion-photon conversion region is filled with a suitable refractive gas. In this 
letter we show that the coupling of axion-photon complex particles to a localized magnetic 
flux generated by a solenoid renders scattering solutions with a cross section which could 
conceivably be measured. 

To see this, let us write the Lagrangian describing the relevant light pseudoscalar coupling 
to the photon. 



Following Ref. {q] (and references therein), we focus on the case where an electromagnetic 
field with propagation along the x and y directions and a strong magnetic field pointing in 
the 2;-direction are present. The magnetic field may have an arbitrary space dependence in 
X and y, but it is assumed to be time independent. 

For small electromagnetic perturbations around the static magnetic background (i.e, the 
axion and the electromagnetic wave), we consider only small quadratic terms in the La- 
grangian for the axion and the electromagnetic fields. By choosing a static magnetic field 
pointing in the z direction and having an arbitrary x and y dependence and specializing 
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to X and y dependent electromagnetic field perturbations and axion fields, the interaction 
between the background magnetic field and the axion and photon fields reduces to 



Ci = -P<j)E, , (2) 

where /3{x,y) = gB{x,y). Choosing the temporal gauge for the electromagnetic field and 
considering only the z-polarization for the electromagnetic waves (since only this polarization 
couples to the axion) we get the following 2+1 dimensional effective Lagrangian 

£2 = ld^Ad>^A + \d^<\>d^<\> - \il<\? + fi<\>d,A , (3) 

where A is the 2;-polarization of the photon, so that = —dtA. 

Without assuming any particular x and y dependence for /3, but insisting that it will be 
static, we see that neglecting the axion mass rua (the validity of this assumption will be 
discussed at the end of this work), we discover a continuous axion photon duality symmetry. 
This is due to a rotational 0(2) symmetry in the axion-photon field space, allowed by the 
axion and photon kinetic terms and by expressing the interaction term, £/, in an 0(2) 
symmetric way by dropping a total time derivative from it: 

Ci = \mdtA - Adt<\>) . (4) 
Defining now the axion-photon complex field, \E', as 

vi/ = _L(0 + ,A) (5) 
and plugging this into the Lagrangian results in 



£ = - -/3(^*5t^ - ^(9t^*) , (6) 



where ^* is the charge conjugation of . From this we obtain the equation of motion for ^ 

d^,d^^ + i^dt^ = . (7) 
We therefore have the magnetic field, or /3/2 (the UiV) charge), coupled to a charge density. 



Introducing the charge conjugation 



10| , that is 
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^ ^ ^* , (8) 

shows that the free part of the action is indeed invariant under ([8]). When acting on the free 
vacuum the A and (p fields give rise to a photon and an axion respectively, but in terms of the 
particles and antiparticles (defined in terms of \E'), we see that a photon is an antisymmetric 
combination of particle and antiparticle and an axion a symmetric combination, since 

= i=(vl/* + VP) and A = -^{^ - . (9) 

Hence, the axion is even under charge conjugation, while the photon is odd. These two 
eigenstates of charge conjugation will propagate without mixing as long as no external 
magnetic field in the perpendicular direction to the eigenstates (i.e axion and photon) spatial 
dependence is applied. The interaction with the external magnetic field is not invariant under 
dH]). In fact, under ([8]) we can see that 

Si -Si , (10) 

where Si = J Cidxdydt. Therefore, these symmetric and antisymmetric combinations, 
corresponding to axion and photon, will not be preserved in the presence of B in the analog 
QED language, since the "analog external electric potential" breaks the symmetry between 
particle and antiparticle and therefore will not keep in time the symmetric or antisymmetric 
combinations. In fact, if the corresponding external electric potential is taken to be a 
repulsive potential for particles, it will be an attractive potential for antiparticles, so the 
symmetry breaking is maximal. 

Even at the classical level these two components suffer opposite forces, thus under the 
influence of an inhomogeneous magnetic field both a photon or an axion will be decomposed 
through scattering into their particle and antiparticle components, each of which is scattered 
in a different direction, since the corresponding electric force is related to the gradient of 
the effective electric potential, i.e., the gradient of the magnetic field, times the U{1) charge 
which is opposite for particles and antiparticles. If we look at the scattering amplitudes for 
particles and antiparticles, we see that they have opposite signs. Calling S the scattering 
amplitude for a particle, the amplitude for an antiparticle is then —5*. Therefore, an axion 
[i.e. the symmetric combination of particle antiparticle (1, 1)] goes under scattering to 
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(1,1) + {S,—S). So the amplitude for axion going into photon (1,-1) is 5". Hence, we 
conclude that the amplitude for axion-photon conversion is equal to the particle scattering 
amplitude. 



For this effect to have meaning, we have to work at least in a 2+1 formalism [ll|. The 1+1 
reduction [3] which allows motion only in a single spatial direction, is unable to produce 
such separation, since in order to separate particle and antiparticle components we need at 
least two dimensions to obtain a final state with particles and antiparticles propagating in 
slightly different directions. 

n 

This is in a way similar to the Stern-Gerlach experiment in atomic physics [12], where 
different spin orientations suffer a different deflection force proportional to the gradient of 
the magnetic field in the direction of the spin. Here, instead of spin we have that the photon 
is a combination of two states with different U{1) charge and each of these components 
will suffer opposite force under the influence of the external inhomogeneous magnetic field. 
Notice also that since particle and antiparticles are distinguishable, there are no interference 
effect between the two processes. 

Therefore an original beam of photons will be decomposed through scattering into two 
different elementary particle and antiparticle components (and also, of course, the photons 
that were not scattered). These two beams are observable, since they both have photon 
components, so the observable consequence of the axion-photon coupling will be the splitting 
of a photon, or axion, beam by a magnetic field of the configuration considered here, whereas 
in the normal Primakoff effect analysis there is no explicit recognition of a splitting. This 
effect is, moreover, of first order in the axion-photon coupling {g), unlike the "light shining 
through a wall phenomena" which depend on the coupling constant squared {g'^)- 



II. FIRST APPROXIMATION: MAGNETIC FIELD OF AN INFINITELY THIN 
SOLENOID 

To apply the results of the previous section to some specific system with magnetic field, 
we write separately the time and space dependence of the axion-photon field as \l/(r, t) = 

e-^'^V(^)- 

As a first model, we are considering an inhomogeneous magnetic field of the form B = 
$(5^(x,y). This kind of a potential can not, of course, be realized in the lab, however, we 
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will show that the results for this, presumably purely theoretical, calculation have physical 
significance in the resonance case, where the scattering becomes isotropic. 

Separating the time and space dependence of \1/ and considering the 6 function potential 
reduces Eq. ([?!) to 

[-V^ + g^E6\x, y)]-^{r) = ^^(r) • (11) 

In terms of momentum space wave functions, 0(/c) = J e^^''^ip{r)d?r, the latter equation 
is now 

P(t){k) + ^$^^(0) = E'^(t){k) , (12) 

from which the solution 

^w = (2.m*-*„)-^^^. (13) 

with kQ = E"^, is obtained. The constant g^Eijj{0) is determined from Eq. (1131) by integration 
over momentum space 

^(0) = 1 - g^Ehi-E^ - ie)V^(O) , (14) 

where 

with and z = —E'^ — ie and A is a cutoff constant that was introduced to regulate the integral 
l2{z) by limiting k. It is straightforward to calculate g^Eip{0) from Eq. ( IT^ 



g^Eip{0) 



1 ^log(AV^) 



g^E Air 



1 ^ log(A/E) ^ 



g^E 2ti 2 



(16) 



To obtain the scattering amplitudes, we write the wave functions in position space 



^(f) = ^ik-r _ g^Eij{Q)Gk{r) , (17) 
where Gfc(r) is Green's function in two dimensions 



(_V2 - k')Gk{r) = 6{r) , 
6 



(18) 



By identifying the scattering amplitude from the asymptotic behavior of the scattering wave 
function 



^P{r) ^ e'^-"''+ i=/(^)e^('^''+"/^) , (20) 



we get for the constant scattering amphtude 

1 Q^E 

since k"^ = E"^. Since there is no dependence on the scattering angle in f{6) the scattering 
from a 6 function is completely isotropic. The total cross-section in 2 dimensions is given 
by atot = Jq^ \fi(^)\'^d6. Hence, by expanding f(6) to first order in g we find that 

= ^ • (22) 

Our primary motivation comes from the QCD inspired axions, with mass up to the ~1 eV 
range. To estimate the magnitude of the total cross-section, we take the value of the coupling 
constant g from the recent result of the CAST collaboration. CAST is searching for axions 
produced in the sun and travelling to earth by trying to detect photons from the conversion 
of axions inside a constant magnetic field, following the coherent inverse Primakoff-effect. 
Along with the Japanese axion helioscope Sumico [13], CAST has set an upper limit on the 
magnitude of the axion-photon coupling constant of < 2.2 x 10~^° GeV^^ for an axion 
mass of nia < 0.4 eV 1^. We choose to use g = 10^^^ GeV^^ throughout this paper. The 
dimensionless magnetic flux is, of course, given by $ = tiBR^, where B is the magnetic field 
strength inside the solenoid and R is the solenoid radius. Lastly, the mean energy of axions 



arriving at the earth from the sun is estimated to be E' = 4.2 x 10^ eV 15 |. 

In order to get the 3D total cross-section (i.e the scattering cross-section) as we multiply 
the 2D cross-section atot by the length of the solenoid L, taking L = 10 cm as an example. 
Multiplying the scattering cross-section by the flux of axions coming from the sun, E = 



3.67 X 10^^ /cm^ ■ sec 15|], we can estimate the number of events per second A^. 

The quantity we are ultimately looking for is the axion-photon conversion probability. 
To obtain this, we calculate the ratio between the number of axions arriving at the solenoid 



to the number of photons produced. The number of axions hitting the solenoid is given by 
multiplying the flux of axions arriving by the geometrical cross section of the solenoid, given 
by ug = DL, where D is the solenoid diameter and L is its length. The number of produced 
photons is found by multiplying the scattering cross section {as = crtot ■ L) times the fiux. 
Thus, the probability is given by 



(^s/(^G 



AD ~ 8 



(23) 



Notice that the dependence on the magnetic field strength is squared. However, the depen- 
dence on the surface magnetic field gradient is "hidden", since it was implied in deriving 
this relation. A few examples for the cross-section, number of events and probability are 
given below in TABLE I. 



B [Tesla] 


D [cm] 


f^L [cm; 




Ns [sec-i] 


Ps = (^s/f^G 


10 


1 


3.08 X 10^ 


15 


0.01 


3.08 X 10~^^ 


10 


10 


3.08 X 10- 


11 


112.98 


3.08 X 10-12 


6 


2 


1.77 X 10- 


14 


0.07 


8.87 X 10-1^ 


6 


20 


1.77 X 10- 


10 


650.78 


8.87 X 10-12 



TABLE I: Total cross-section, number of events and axion-photon conversion probablity for 
different choices of the magnetic field strength {B) and the solenoid diameter {D) and for 
g = 10-1" GeV-i. We have used rationalized natural units to convert the magnetic field units 
from Tesla to eV^, where the conversion is 1 T = 195 eV^ (please see appendix A in 161] for more 
details). 



III. FINITE SIZED SOLENOIDAL GENERATED POTENTIALS 
A. Gaussian Distributed Magnetic Field 

We wish to obtain eventually measurable quantities which can be incorporated in a labora- 
tory experiment, thus we have to consider a more realistic function to describe the magnetic 
field generated by the solenoid. As a first model, we choose to describe the inhomogeneous 
magnetic field by a Gaussian distribution around the solenoid's axis. 
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^2 



B{r) = Boe^z . (24) 

Introducing again Green's function in the x , y plane, we write the wave function in 
position space 

m = ^/r-ee(r) + j G{r- p) Q E B {p),p {p) dP p , (25) 

where il^free = e*'^'^ is the solution of the free field equation. To first Born approximation, 
noting that 



and using again the asymptotic approximation of Green's function in 2 dimensions (see Eq. 
[T9l) we arrive at 

V^(r) = e'^^-^ + -^== / gEBip)e^^^-^^d'p , (27) 

where q = k — k^. To evaluate the integral, B{q) = J B{p)e'^^'^cfp, we write p-q = qp cos (j) 
and get 

/•OO 2 /'2'K l*00 2 

Bo / eiJ^pdp / d^e'^P'"""^ = 2tiBq / e^Jo(gp)prfp = TiBoR^e--^^^^^' . (28) 
Jo Jo Jo 

Hence, the wave function becomes 



^ ^^k.r ^ ^/^9Bo^^^-liR,f^^ikr+n/4) (29) 



By defining, as before. 



V^(f) ^ e'^-"> ^/(^)e^(fc'-+-/4) ^ (30) 



we find for the scattering amplitude 

f{e) = y^)gBoR'E'/'e-'^^^^^' , (31) 
where the explicit dependence of q on the angle is given by 



Hence, The total 2D cross-section is given by 



(32) 



\f{d)\^dd 



TT 



Jo ^ 



(33) 



where Io{x) = Joiix) is the modified Bessel function. The argument of this function (i.e 
{RkY) is very large (1 eV x 1 cm 10^) so we can use the asymptotic from of the modified 
Bessel function 



1 + 



;i-2n)(l + 2n) 



y/2TXX V ' Sx 
Keeping only the first order term gives the result 



+ ... 



(34) 



Gauss 
tot 



IT 



3/2 



32^ 



(35) 



Again, we find the axion-photon conversion probability P = as/o'c to be 

7,3/2 



Gauss 



8V2 



g'BlR^ . 



(36) 



a result which is about two times larger than the ID case 



171] (when taking the linear 



dimension associated with the extent of the magnetic field as the solenoid's radius). 



B [Tesla] 


D [cm] 


agr^^ [cm] 


Ncauss [sec 




PGauss 




10 


1 


1.17 X 10"23 


4.29 X 10" 


11 


1.17 X 10" 


-23 


10 


10 


1.17 X 10^20 


4.29 X 10" 


-8 


1.17 X 10" 


-21 


6 


2 


3.38 X 10"^^ 


1.24 X 10" 


10 


1.69 X 10" 


-23 


6 


20 


3.38 X 10^20 


1.24 X 10" 


-7 


1.69 X 10" 


-21 



TABLE II: Total 2D cross-section, number of events and the axion-photon conversion probablity 
for different choices of the magnetic field strength (B) and the solenoid diameter (D) for the finite 
sized solenoid with Gaussian distributed magnetic field case. Again, we use g = 10"^*^ GeV~^ 
and rationaUzed natural units to convert the magnetic field units from Tesla to eV^, where the 
conversion is 1 T = 195 eV^. 
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B. Solenoidal Generated Potential - Square Well Approximation 



Now we turn to consider the magnetic field generated by an ideal solenoidal current which 
is described by a step function realizing a uniform magnetic field pointing in the z direction 
and constrained to a cylindrical region around the origin 

Bnz , r < R , 
B{r) = { (37) 

, r>R. 

Repeating the same manipulation as in equations fl2Sl) to fl52|) and using the Fourier trans- 
formation of the step function 

Bo [ pdp [ d^e'''P'°''^ = 2nBo [ pdpJo{qp) = ^^^^Ji(gi?) , (38) 
Jo Jo Jo 1 

we find that the scattering amplitude is now given by 

m ^ f^?^U,n) . (39) 

where the explicit dependence of q on the angle is given by Eq. (132|) . 

Before evaluating the integral for the total cross-section, let us write the total cross section 
for the square well case in terms of the delta function cross-section, calculated in section [III 



a 



well 
tot. 



IT 



—g'BWE 



2n 



MqR) 



qR 



de 



27r 



Ji{qR) 



qR 



dO 



al,-I{ER) 



(40) 



where I{ER) = /J 



2lT 



■h{qR) 



qR 



dO is a dimensionless quantity which is a function of the mul- 
tiplication E ■ R. Using the relation P = atot/D (where we use the same notations as in 
section [IT]) , the proportionality constant connects also the conversion probabilities for the 5 
function and square well cases 



well 



Ps-I 



g^B^D^EI{ER) . 



TT 32^ 

Denoting ER = kR by t], the integral can be analytically solved with the solution 



(41) 



TT 



I{v) = ^2Fs{{il}; {1,2,3} -,-4^' 



(42) 
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where 2-^3 is an hypergeometric function. 

To analyze this solution we expand the hypergeometric function to a series. Then, for 
small rj, I{r]) is converging toward the constant value 7r/2, thus giving the equality cr^f' = 
crfgj . This result is expected since considering only small r] values is equivalent to considering 
isotropic scattering because ^ 1 means that ER ^ 1. Hence, the wavelength of \E' is 
very large compared to the length scale of the potential. Therefore, this approximation 
corresponds to 6 function limit of the step function, which, in turn, means that we consider 
isotropic scattering. 

This conclusion can also be deduced from the following viewpoint regarding the scattering 
angle: The integrand of / is becoming extremely oscillatory as its argument (i.e. qR) is 
bigger and therefore for a reasonable scale of ER 10^) we have a highly oscillatory 
integrand which is also decaying very fast as a function of 6 (since the momentum transfer q 
is a function of the scattering angle). Thus, the biggest contribution will come from smaller 
angles. In fact, demanding that the integrand will be of order one is equivalent to considering 
scattering ang les that satisfy 6 < 1/ER ^ 10"^ Then, using the asymptotic form of 
the Bessel function for small arguments we have I ^ 7^/2 which simply gives the solution 
^Tof" — ^tot — T^"^ g"^ B"^ E / A. Considering only small angles is equivalent to demanding that 
the argument of the Bessel function will satisfy the condition ER ■ sin(^/2) ^ 1. Without 
limiting the range of the scattering angle, this of course means that 77 = ER <^ 1, the 
condition which coincides with the small rj expansion of I{ri). 

On the other end, we have the expansion for large rj. This reveals the fact that the 
integral approaches the limit / — > = 3^;^ very fast. For example, for r] = 10 we already 
have -^/I{ri = 10) = 0.997. A plot of I{ri) and its limit is shown in Fig. [1] Putting 
this limit into Eq. fj4T]) gives the result 

Pn^eii = Ig'B^D-' = 2|PiB , (43) 

1 1 

where Pid = ^g^B'^R^ is the ID conversion probability [15]. Thus, the scattering from a 
step function potential enhances the probability of the ID case by a factor of 2.67. 

Since the generalized hypergeometric function is difficult to analytically work with for 
large arguments (large 1] values), we have also calculated the total 2D cross-section numer- 
ically to verify our results for the entire spectrum of 77. In order to evaluate the integral / 
we have used the 'MATLAB' program, running the new 'quadgk' function which is using 

12 




FIG. 1: The solution of the integral /(??), defined by Eq. HOl is an hyper geometric function I(rj) = 
f 2-p3({|5 §}; {1, 2, 3} ; — 4ry^). This figure shows a plot of /(r/) as a function of the multiplication 
ER = 7] in the solid line. The dashed line represents the fast approached limit of /(??), which is 
given by At 77 = 10 both lines are close enough so that the ratio 55;; /-^(^ = 10) already equals 



0.997. 



the Gauss-Kronrod quadrature and is efficient specifically for oscillatory integrands. The 
results of the the analytical and numerical calculations match to very high precision and 
both results are practically the same. In fact, when considering solar axions (i.e. ER = rj 
is of order 10®), the numerical calculation gave the result Pweii = 2.67Pi£) as well. A few 
examples for the cross-section, number of events and conversion probability for solar axions 
are given below in TABLE III. 

When comparing the cross-sections of the Gaussian distributed magnetic field to the 
step function, one expects the step function cross-section to be bigger than a cross-section 
generated by a smooth function, in agreement with similar studies done in the context of 
nuclear physics models, where it has been shown that a step function potential gives a bigger 
scattering cross-section than a smooth potential like, for example, the Woods-Saxon Diffuse 
potential 18|. Our results qualitatively agree with Woods and Saxon, as can be seen by 
comparing Tables II and III. 
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B [Teslal 


D [cml 

L J 


aZ^^ [cm] 

ZOZ L J 


N,„f,ij [sec" 


''] 






10 


1 


1.58 X 10^23 


5.80 X 10" 


11 


1.58 X 10^ 


'23 


10 


10 


1.58 X 10-20 


5.80 X 10- 


-8 


1.58 X 10" 


-21 


6 


2 


4.56 X 10-23 


1.67 X 10- 


10 


2.28 X 10" 


-23 


6 


20 


4.56 X 10-20 


1.67 X 10" 


-7 


2.28 X 10" 


-21 



TABLE III: Total 2D cross-section, number of events and the axion-photon conversion probablity 
for different choices of the magnetic field strength (B) and the solenoid diameter (D) for the finite 
sized ideal solenoid case. We use g = 10"^*^ GeV-^ and rationalized natural units to convert the 
magnetic field units from Tesla to eV^, where the conversion is 1 T = 195 eV^ 



IV. RESONANT SCATTERING FOR S ~ 

So far in this report, we have consider the axion field as a massless field in order to get 
the U{1) symmetry between axions and photons. In fact, this symmetry holds up whenever 
the axion mass is equal to the (effective) photon mass inside a medium. For example, in 
axion helioscope experiments photons acquire an effective mass if one fills the conversion 
region with a suitable refractive gas. 

Of course, if we had recalculated our results with massive axions and "massive" photons 
(of equal mass to that of the axion), our conclusions will have to be modified. The term 
that has to be taken under consideration is an 1/ a/ (-E^ — m^)!/^ term which comes from the 
Green's function and will replace the current 1/\/E in the scattering amplitude. Thus, in 
the rria ~ mphoton 7^ case, the total two dimensional cross-section (for the 6 function case) 
would have the following energy dependence 

^tot = . „ = , (44) 

and we have a resonance when E = m, which has, in a sense, a similar behavior to the ID 



problem analyzed by Adler et. al. [21[ (notice that Adler et. al. consider the conversion 
between a massive axion and a massless photon), where of course the resonance here appears 
at ruaxion ~ rriphoton- In Eq. (jSj) the relation rriaxion ~ rriphoton is assumed from the beginning 
and we see that the additional resonance appears as ~ m. For an axion rest mass below 
~ 1 eV, this can have practical consequences, for example, in laser generated axions (e.g in 
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'shining through the wall' experiments) when one can control the energy of the axion beam. 

We can see here that the ID treatment of this process can not be justified since in the 
limit of zero momentum the scattering amplitude and the differential scattering cross-section 
become isotropic (i.e equal for all angles) and it is impossible to consider only one direction 
in the scattering. In fact, in the limit of exactly zero momentum (assuming there is a tunable 
laser capable of very fine accuracy to obtain E very close to m) the amplitude of a finite 
potential becomes of the form to Eq. (1441) . This is since taking the limit of zero momentum 
implies zero momentum transfer (from Eq. |32|) which means to consider only zero modes 
in the Fourier transform of the magnetic field. Hence, the cross-section of a finite potential 
becomes of the form of the modied delta function potential. It is an experimental question 
whether such a fine tuning is possible with an existing laser, if the axion has a mass of the 
order of eV. 

Achieving a resonance requires a material which has a zero index of refraction. The real 
part of the refractive index is given by 

nn{uj) = l + K ""'-^ ^ , (45) 

where K = Ne^f with N being the number density of atoms, e is the electron charge and 
/ transition oscillator strength, is the transition frequency and 7 represents dissipative 
interactions [19] . Equating the latter to zero requires the condition > 47^. A negative 
and zero refractive indices are indeed possible as was experimentally observed by Shelby et 



al. 20| . Let us hope that one day it will be possible to implement this in an axion detection 



lab experiment. 



V. SUMMARY, DISCUSSION AND CONCLUSIONS 

In this paper we have studied the first examples of scattering which is not one dimensional 
and we have obtained enhanced probabilities. This effect is further increaed in the case of 
resonant scattering that appears when E = m and corresponds to isotropic scattering (as in 
the 5 function scattering). One should notice that allowing for two dimensional scattering 
is the same as allowing the possibility of axion-photon splitting which does not make sense 
in ID scattering. We have studied here merely magnetic fields with a cylindrical structure. 
Further generalizations should include the scattering from, for example, a quadrupole mag- 
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netic field, which is more complicated than the cylindrical symmetric case we have studied 
here but, on the other hand, is quite accessible as a possible experimental setup. 

In the ID case the conversion probability is Pio = g'^B'^P /A 17|, where / is the linear 
dimension associated with the extent of the magnetic field {Pm = Pcauss./ i'^'^'^^'^))- The 
comparison between the ID scattering, finite sized potential scattering (for example, given 
in the table is the step-function potential) and isotropic scattering {6 function case) is given 
below in TABLE III. Looking at this table, we see that the effect of considering 2D scattering 
instead of ID scattering increases the probability for an axion to be converted into a photon 
for axions coming from the sun. In fact, we see that Pweii = 2.67 x Pid (as was shown in 
Sec. UTTBD . 



B [Tesla] 


D [cm] 


PlD 




Ps 




Pwell 




10 


1 


5.94 X 10- 


-24 


3.08 X 10" 


-15 


1.58 X 10" 


-23 


10 


10 


5.94 X 10" 


'22 


3.08 X 10" 


-12 


1.58 X 10" 


-21 


6 


2 


8.56 X 10" 


-24 


8.87 X 10" 


-15 


2.28 X 10" 


-23 


6 


20 


8.56 X 10" 


-22 


8.87 X 10" 


-12 


2.28 X 10" 


-21 



TABLE IV: A comparison of the axion-photon conversion probability of the ID axion helioscope 
case (-Pid), the 2D delta function case (P^, calculated in section [Tl]) and the 2D step function {P2D, 
calculated in section ITlI Bp . We use, as in the rest of the paper, g = 10~^^ GeV~^ and rationalized 
natural units to convert the magnetic field units from Tesla to eV^, where the conversion is 1 T = 
195 eV^ 



When considering scattering from a finite sized potential (Gaussian and step function 
potentials) the enhancement of the conversion probability compared to the ID case still gives 
probabilities in the same order of magnitude. This is due to the fact that the wavelength 
(l/E) of the \l/ wave function is much smaller than the length scale of the potential {R), 
which essentially results in a quasi-lD behavior of the system. When the wavelength will 
be smaller, or even comparable to the length scale of the potential we see that we get bigger 
enhancement since in this case the scattering becomes more and more isotropic and we 
essentially obtain S function scattering. 

The wavelength is determined by the momentum of the particles (from the de Broglie 
relation). Hence, the smaller the momentum the bigger the wavelength. For the massive 
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case, the momentum approaches zero when the energy of the particles is of the order of the 
particle's mass. This situation, where the wavelength of the particles is much larger than 
any other length scale in the problem, is realized in the resonant scattering case, discussed 
in Sec. IIVI There we have shown explicitly that this limit gives an isotropic scattering for a 
finite potential and thus, conversion probabilities of the order of the S function case (shown 
in TABLE I). 

The cross-section in the resonance case was calculated at tree level. This gives a singu- 
larity of the cross-section at E = m. However, in practice a resonance effect should have 
a certain width and this, of course, should also be the case for the resonance case found 
here. We notice also that the resonance behavior comes together with a breakdown of the 
ID treatment of axion-photon conversion and also that a finite width can be originated from 
absorption effects. All these problems will be addressed in a future publication. 

Our results might also be applicable for the solar scenario as well. In the sun, magnetic 
flux tubes can play the role of a solenoidal potential while the energy spectrum of photons is 
continuous. Thus, we expect to have both isotropic (resonance) and anisotropic scattering. 
These magnetic flux tubes are enormous regions of constant magnetic flux with length scale 
of the order of about 10^ km in diameter and 10'^ km in length. If we trust our numerical 
results to work in these scales as well, the conversion probability (which of course relates 
to at-least 2D scattering) will be greatly enhanced. For example, taking a flux tube with 
magnetic field of -B = 0.2 T and diameter of 100 km, we get a conversion probability 
of Psun = 6.34 X 10-13 g ^ iQ-w GeV^i), which is even larger than the isotropic 

scattering from a laboratory fictitious infinitely thin solenoid. 

This result may be related directly to some basic solar physics problems and, in particular, 
the coronal heating problem in the sun. The sun's outer layer, the solar corona, is much 
hotter than the surfaces below it, the chromosphere and the optical surface of the sun (the 
photosphere). Within a few hundred kilometers, the temperature in the corona rises to be 
about 500 times that of the underlying chromosphere, instead of continuing to fall to the 
temperature of empty space (2.7 K). While the energy flux of extreme ultraviolet photons 
and X-rays from the higher layers of the sun is some five orders of magnitude less than 
the energy flux from the photosphere, it is still surprisingly high and inconsistent with the 
spectrum from a black body with the temperature of the photosphere. Thus, some exotic 
physics must be at work out there. We would like to point out here that even the rather 
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modest probability enhancements derived in this work, of as much as a factor of 2.67, might 
still provide a potential explanation we are looking for, as for why the solar X-ray activity 
correlates preferentially not only to magnetic fields but even more so in places with magnetic 
field gradient (near the inversion line of two oppositely directed magnetic field regions) 22 1. 

Moreover, as was shown earlier in this work, when the axion mass and the energy corre- 
sponding to the plasma frequency are equal the conversion probability features a resonance 
and increase sharply. Since in the restless magnetic sun the magnetic fields and plasma 
densities are continuously changing, this resonance crossing is quite probable and can result 
in an otherwise unexpected photon excess or deficit. Hence, the work presented here might 
suggest, in its astrophysical extension, a possible solution to some basic problems in solar 
physics. These issues will continue to be studied by us in the future. 
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